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AN APPLICATION OF THE INVERSE SCATTERING 
PROBLEM FOR THE DISCRETE DIRAC OPERATOR 

RZA I. ALESKEROV 


Abstract. The method of the inverse scattering problem is used to 
solve the Cauchy problem for the Langmuir lattice with initial data of 
step type. Formulas are obtained for transforming the scattering data 
with respect to the time, making it possible to obtain a solution of the 
problem for arbitrary moment with the aid of linear integral equations 
of scattering theory. The existence of a solution is proved. 


1. Introduction 


It is known that the method of the inverse scattering problem allows a detailed 
study of the Cauchy problem for some nonlinear evolution equations (see [1], 
[18],[19] and the literature therein). In works [2]-[7], [9]-[11], [17],[20] the Cauchy 
problem was studied for nonlinear difference equations with different initial data. 
In solving these problems, inverse spectral problems are used for the discrete 
Sturm-Liouville operator. On the other hand, the inverse scattering problem 
for the discrete Dirac operator was studied in [12]-[15]. However, as far as we 
know, applications of the latter problem to nonlinear equations have not been 
considered. In this paper, we apply the inverse scattering problem for some 
discrete Dirac operator to the Cauchy problem of the Langmuir lattice with a 
step-like initial condition. 

It should be noted that the Langmuir lattice has important applications in 
plasma physics and in zoology. The Cauchy problems for this lattice in various 
classes of initial conditions were studied in [2],[11], [17],[18],[20]. 

For a sequence of positive functions c n = c n ( t ), c n £ [0, oo) we consider 

the Langmuir lattice 


C-n 


Cn (Cn— 1 Cn+l) , • 


A. 

dt 


n £ Z. 


( 1 . 1 ) 


For (1.1) we pose the following Cauchy problem with the initial condition 


c n (0) = c°,ra £ Z, 


( 1 . 2 ) 
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satisfying requirement 

— 1 OO 

^2 M (l c n (0) - 1|} + ^ |n| {|c n (0) - A|} < oo, (1.3) 

ii=—oc n= 1 

where 0 < A < 1. 

We will seek a solution to problem (1.1) and (1.2), such that 


-i 

E 9 + N) 

n =—oo 


C n (t) - 1| + ^2 (1 + l n D \ C n (*) “ ^1 

n=0 


< OO 

C[0,T] 


(1.4) 


for all T > 0 . 

In this paper, based on the Gelfand-Levitan-Marchenko formalism, the problem 
(1.1)-(1.2) is studied. Explicit formulas are obtained that describe the evolution 
of scattering data. The unique solvability of problem (1.1)-(1.2) in the class of 
functions satisfying (1.4) is proved. We note that the asymptotic behavior of the 
so-called rapidly decreasing solution of problem (1.1)-(1.2) using the Riemann 
problem was studied in [20] . 


2. THE UNIQUE SOLVABILITY OF PROBLEM (1.1) - (1.2) 

Theorem 2.1. Problem (1.1)-(1.2) has a unique solution in the class of functions 
satisfying (1-4)> if the condition (1.3) is satisfied. 


Proof. Introduce the Banach space E of sequences x = (x n )^_ 0O with the norm 

OO 

IMIe = (1 + l n l) \ x n\- Then, the set C([0, T]; E) of functions continuous on 

n =—oo 

[0, T] with the values in E is a Banach space with the norm ||x (i)llc([o t]-e) = 
||x (t)\\ E (see, e.g.,[16] ). Assume 


X n (t) 


c n (t) - A, n> 0, 
c n (t) - l,n < 0. 


Then the problem (1.1)-(1.2) is equivalent to the problem 

x n = x n (x n -i - x n+ i) + A (x n _i - x n+1 ) ,n> 1, 
x 0 = x 0 ( x-i - xi) + A ( x-i - xi) + (1 - A) xq + A (1 - A ), 
X-l = X-l (x -2 - x 0 ) + (x _2 - X-l) + (1 — A) x 0 + (1 — A ), 
x n — x n {Xn— 1 ®n+l) T {.Xn— 1 ®n+l) > ^ ^ 1- 


( 2 . 1 ) 


( 2 . 2 ) 


x n (0) = x°,n G Z 

where 

o = / c° - A, n > 0, 
n \ c°-l,n<0. 

Rewrite problem (2.2) and (2.3) in the form of an operator equation: 


T 


x ( t ) = x (0) + 


T (x (t)) dr 


(2.3) 


o 
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where x(t) = (x n (t))’^ = _ oc and 4' is the operator generated by the right-hand side 
of system (2.2). It is obvious that, for each T > 0, the operator is a continuously 
differentiable mapping of the space C([0,T]; E) into itself. Therefore, we can 
apply to the last equation the contraction mapping principle. As a result, we 
find that, in a certain interval [0, 5], problem (2.2) and (2.3) has a unique solution 
x(t) = (x n (£))“L_ 0O and with a finite norm \\x (i) IIo([o 5]-E) • Let us P rove that 
this solution can be continued to the finite interval [0, T\. Assume the contrary. 
Then, there is a point to < T such that problem (2.2)-(2.3) has a solution x(t) = 
(xn(t))n l-oo continuous in the interval (0,to) 5 but lim ||x(f)|| E = +oo. 

t-*to 

Using relations (2.1), we obtain 


(2.4) 


x n ( t ) = x n (0) + / [x n (r) (x n -i (r) - x n+ i (r)) + 

o 

+ A (x n -i (r) - x n+ i (r))] dr, n > 1, 

t 

Xo (t) = xq (0) + A( 1 - A)t + f [x 0 (r) (x_i (r) - xi (r)) | + 

o 

+ A (x-i (r) - xi (r)) + (1 - A) x 0 (r)] dr, 

t 

X—i (t) = X-! (0) + (1 - A)t + / [x_i (r) (x -2 (t) - x 0 (r)) + 

0 

+ (x —2 (t) - x_i (r)) + (1 - A) x_i (r)] dr, 

t 

x n (■ t ) = x n (0) + / [x n (r) (x n _i (r) - x n+ i (r)) + 
o 

+A (x n _i (r) - x n+1 (r))] dr, n < -2. 

It is known (see ) that problem (1.1) and (1.2) has a unique solution bounded 
uniformly with respect to t in each finite interval, i.e., \c n (t)\ < C,n G Z,t G 
[0, T}. It follows that the problem (2.2)-(2.3) also has a unique solution x n (t) 
bounded uniformly with respect to t in each finite interval: |x n (£)| < M,n G 
Z,t G [0, T] . Using the last fact, from relations (2.4) we have 

t. 


E < II* 


E + (5 M + 1) T + (2 M + 2) / ||x (r)\\ E dT, 0 < t < to. 


o 


By the Gronwall lemma, we find from the last inequality that 

]l* m E < 01* m\ E + (5 M + 1) T] exp ((2 M + 2) T ), 

which contradicts our assumption. Hence, problem (2.2)-(2.3) on the interval 
[0, T] has a unique solution with a finite norm ||x (i)llc?([o t\-e)- ^he theorem is 
proved. □ 


3. EVOLUTION OF THE SCATTERING DATA 

Consider problem (1.1) - (1.2). Putting in equation (1.1) 

C2n — ®l,n) C2n— 1 — ^ 0,(l2,n < 0, Tl G Z , 


(3.1) 


ai,n = \dl ,n {4,n-l ~ a 2 ,n) = 

^2 ,n = 2 ^ 2 , in (^l,n ^l,n+l) ? ^ ^ 


we obtain 


(3.2) 
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a l,n (0) — 0, l n , 0-2,n (0) — 0,2,ni (3-3) 

where ai )n ,a 2)Tl satisfy the condition 

— 1 OO 

y, {|n| (|ai jn - 1| + |a 2 ,n + 1|)} + ^ (M (\oi, n - A\ + |o 2 , n + A|)} < oo. 


n=l 


According to the results of the previous paragraph, problem (3.2)-(3.3) has a 
solution ai : „, = ap n (f), 02 , n = 02 , n (t ) satisfying the condition 


E (1 + |n|) {|ai,rc ( t ) — 1| + 1 02, n (t) + 11} + 

n =—oo 
00 

+ X/ (1 + |n|) {|fli,n (t) — A| + \a 2 ,n (t) + A|} 

n =0 


(3.4) 


< 00 


C[0,T] 


for all T > 0. 

Now we consider a system of equations 


0\,nV2,n 02,nV2,n+l — 

Ol,n—iyi,n—l 4 “ 02,nUl,n — ^y2,in'0 £ 


(3.5) 


in which the coefficients ai )2 ,a 2;n depend on t and are a solution to problem 
(3.2)-(3.3) in class of functions satisfying (3.4). For definiteness, we assume that 
A < 1 . For solving problem (3.2)-(3.3), we use the scattering theory developed 
in [13],[14] for the system of equations (3.5). 

Let Tj be the complex A- plane with the cut along segment [—A 2- - 7 , A 2- - 7 '] and 
denote its boundary j = 1,2. Consider the function 


Zj = Zj (A) = - 


A 2 - 2 A 2{ - 2 ~ti A 


2 AW~i) + 2A^7^-dA^ 2 J), 


where the regular branch of the radical is selected from the condition \J\ 2 — 4A 2 ( 2 A 
0 for A > 2 A 2 ~i,j = 1,2. As follows from [13],[14] the system of equations has 

solutions | f+ n (A,t) j and j f~ n (A,f) j, j = 1,2, representable in the form 


ftn ( A > 0 ( A > 0 = a j ( n > f ) (yy ) 2 J z i y 1 + K t ( n > m > t ) z Tj > 

fj, n ( A > 0 = a j (”. t) (yy ) z 2 n \ 1 + E K i (n, r >h t)zl m ] , n G Z, 

' ' \ m <—1 / 


where the quantities (n, t), (n, m,t) ,j = 1,2 satisfy the relations 

y (n, f) = 1 + o (1) n —>- ± 00 , j = 1, 2, 

K± ( n , to, t) = O (cr ± (n + [??] + , n + m —±00 




(3.6) 


here a 1 * 1 (n) = E 1 


. Ml 

0\,m A 2 


+ 


« 2 ,»n + A 2 | and [x] is the integer part 


zbm>zbn 

of x. Substituting the representation (3.6) of the Jost solution in the system of 
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equations (3.5), we have 


a l,n 

TIT 


af (n+l,t) \ 
af(n,t) ) 


±1 


a 2, 5 

i±i 


_ ( 

f(n,t)) 


A 

a l’ n A i±i i±i = ± (Ki (n+ ^,±1 ,t) -Kf (n + ^,±l,t)) , 

= ± [Kf (n - 1 1 1 , ±1, t) - Kf (n + ±fl, ±1, t)) , n G Z 
These solutions are coupled on the continuous spectrum by the equalities 


(3.7) 


f~ n (A, t) = a, (A, t) /+ (A,t) + h (A, t) /+ n (A, t) , A g cdd, A 2 jk 4A 2 , 
/,+ (A, t) = G 2 (A, t) f- n (A, t) + b 2 (A, t) f- n (A, t), A g dV 2l A 2 + 4, 


(3.8) 


where the coefficients aj (A, f), bj (A, t) ,j = 1, 2 are determined by the formulas 


«i 


a t) = v)] 

^ ’ ’ A 2 (zi—z 1 J ) 


bi (A,i) = 


AVE 


a 2 (A,i) = 




( 2 : 2 -z 2 ) 


MA. t) = w ^ (v) ^ (v > 


( z 2 Z2 ) 

The coefficients ai(A,t), a 2 (A, t) are limiting values of functions analytic in the 
plane T 2 , wherein they have a finite number of simple zeros A& = ±/i*. = 
±/q,. (t), Hk > 0, k = 1,..., IV. Assume 

12 ' 


(m+ (t)) 2 = { | /£» (±A*k> *) + fin t) }, (A, i) = 

nez ' ' 


K (*)) 2 = X] 1 |/l,»(=■=»>*) + An (=*=». *) I !><• (A, *) = 


£ 

n£Z 


h CM) 

ai (A,f) 

b 2 (A, f) 
a 2 (A,f) 


where we took into account that (see [14] ), to symmetric eigenvalues T/j/,, corre¬ 
spond equal normalization coefficients rnf (t) and m ^ (f). 

We call the set {r + (A, t) , A G <9Ti; Hk (t ); mf (t) ,k = 1, 2,..., IV } and 
{r“ (A, f),A £ <9T 2 ; Hk (t ); m]T (t) ,k = 1, 2,..., IV} the right and left scattering data 
for the system of equations (3.5) respectively. It is appropriate to emphasize 
that the right scattering data is uniquely determined by the left scattering data. 
The inverse scattering problem for the system (3.5) is recovering the coefficients 
Gi )2 , CL 2 ,n by the left scattering data. From the left scattering data we can recover 
the coefficients ai i2 ,a 2;n in (3.5). Let 

N „ . / -1 a ,2(2 -j) 


F 3 M) = Z K W) 2 (^X 2 ) 


1 _T_ f Ar A4) ( z 2 - 1 

2-rri J < 1 


dr 2 


O^ 1 -**) V A 


2(2—J) 


A = ±i2 fe 


(3.9) 


l dA. 


The kernels K 1 (n, rn, t) , iv 2 (n, m, t) are obtained from the main equations of 
the Gelfand-Levitan-Marchenko type 


Kj (n, m, t ) + Fj (2 n + m, t) + 
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+ ^2 K 7 ( n ’ r , 0 ( 2 n + m + r,t) = 0 ,m< -1 ,j = 1 , 2 , (3.10) 

r< —1 

which have unique solutions for arbitrary n. Following this, we determine the 
coefficients 01 , 2 , 02 , n from either one of the formulas (3.7). 

Theorem 3.1. Consider the system of equations (3.5) with coefficients ai, n = 
ai,n(0, a 2,n = a 2 ,n(f) ■ If these coefficients 01,2,02,71 make up a solution of equa¬ 
tion (3.2) in the class of functions satisfying (3.f), then the evolution of the left 
scattering data is described by the following formulae: 

r~ (A, t) = r~ (A, 0) exp { ( zf 1 — z 2 ) t} 

Vk (t) = n k (0) = ia k , (3.11) 

( m k( t )) 2 = {rnf( 0)) 2 exp { [zf 1 (qi k ) - z 2 (//*,)) t} , k = 1,..., N 

Proof. We consider the operator L = L(t) generated by the right-hand side of 
the system of equations (3.5). We also introduce the operator A = Aft), setting 

(^2/)l.n = l®2,nj/l,n— 1 3 ®l,7i02,n+l2/l,n+l 1 

(■^■2/) 2. n — 2®l> n— l®2,n— lU 2 ,n— 1 2 ®2,n2/2,n+l • 

It is easy to verify that the operators L and A form a Lax pair, and therefore the 
system of equations (3.2) is equivalent to the operator equation 

L = AL — LA. 

By virtue of the last relation, the operator B = 7 — A takes solutions of the 
system equation with parameter t to solutions of the same system. Using (3.6), 
we can show by standard methods (see ) that, as n -A 00 

A3 / A r _ A \ 

[Bf + (A ,t)] jn = — (zf 1 - zi) y—J -J z i + o(l) ,n -> 00 

A solution with such asymptotic behavior is unique; hence, 

[Bf + (A, t)] jn = (zf 1 - Zl ) ff n (A, t),j = 1 ,2 

On the other hand, from (3.8) we obtain 

Bf + (A, t)] jn = [d 2 (A, t) + \ (zf 1 - z 2 ) a 2 (A, t)] f~ n (A, t)+ 
b 2 (A, t)-\ (zf 1 - z 2 ) b 2 (A, t ) f~ n (A, t),j = 1, 2 . 

Comparing the last equalities with (3.8), we obtain 

a 2 (A, t) + 3 (z ^ 1 - z 2 ) a 2 (A, i) = (zf 1 - zi) a 2 (A, f), 
b 2 (A, t) - 3 (z ^ 1 - z 2 ) 62 (A, i) = 4“ (%* - 21 ) b 2 (A, t). 

Whence it follows that 

a 2 (A, t) = a 2 (A, 0) exp j 4^ (zf 1 - zi) t - | (z ^ 1 - z 2 ) t| , 

6 2 (A, f) = fe 2 (A, 0) exp 14r (zf 1 - Zi) t + 3 (z ^ 1 - z 2 ) tj • 

As follows from last formulas, the zeros of the function 02 (A, t) do not depend on 
t and the first two relations in (3.7) are valid. 

We now find the law for transforming the normalization mf (t). 
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Let ipj t n (^k,t) be a normalized eigenfunction of the operator L corresponding 
to the eigenvalue /i^: 

OO 

£ {Vl,n(Vkit) + V2,n(Vk,t)} = 1, 

n =—oo 

Pj,n (pk, t)~c (t) ( ^ )-1 ) z 2 n (Mfc), n -> -oo 


Since 


fj,n 


z 2 (/hfc) - 1 




2 -j 


z 2 71 (jUfc), n —> -oo 


it follows that 

fj,n (Pki t ) = [l 1 ki t) , 

and consequently 

W )” 2 = S { A" n + / 2 7 n | = c 2 (t). 

™ez ^ J 


It is easy to verify that the normalized eigenfunctions <*p hn (nk,t) satisfy the equa¬ 
tion 

Pj,n (pk, t) - (Aip (n k , t)) jn = 0 ,j = 1 , 2 . 

Therefore, by virtue of (3.12), we have 

c (t) - ^ (z 2 x (ji k ) - z 2 (nk)) c(t) = 0 

whence we find that 

c (t) — c (0) exp | * (z 2 1 (/x fe ) - £2 (/ifc)) * j 

( m fc (*)) ^ ( 0 )) _2 exp {(/ifc) - z 2 (/ik)) tj 

This completes the proof of the theorem. □ 


The theorems proved above give an algorithm for finding a solution to prob¬ 
lem (1.1) - (1.2). For constructing the solution of problem (1.1)-(1.2) by the 
initial data ayn(0), j = 1, 2 , we calculate the scattering data {r~ (A, t ), A G <9r 2 ; 
Hk ( t ) \m ]7 (t) ,k = 1,2, ...,1V}, the of the system of equations (3.5) for f = 0 . 
Then the set {r~ (A, t), A £ <9r 2 ; Hk (t ); m )7 (t) ,k = 1, 2,..., IV} may be found 
by means of (3.11). Construct the function Fj~(n,t) by the formula (3.9). Find 
the solution of the equation (3.10). Calculate c n = c n (t) by the formula (3.1), 
(3.7). 
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